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Oh' Abstract. We have carried out a survey of the Andromeda galaxy for unresolved microlensing (pixel lensing) . We 

present a subset of four short timescale, high signal-to-noise microlensing candidates found by imposing severe 
' selection criteria: the source flux variation exceeds the flux of an J? = 21 magnitude star and the full width at half 

^ I maximum timescale is less than 25 days. Remarkably, in three out of four cases, we have been able to measure or 

. strongly constrain the Einstein crossing time of the event. One event, which lies projected on the M31 bulge, is 

almost certainly due to a stellar lens in the bulge of M31. The other three candidates can be explained either by 
stars in M31 and M32 or by MACHOs. 
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Key words. Galaxy: halo - M31: halo - lensing - dark stars. EROS IjLasserre et al. 2000|l obtained upper limits 
matter that are consistent with the MACHO results. 

Since M31 is 15 times more distant than the Magellanic 
Clouds, the stars are about 200 times fainter and more 
1 Introduction densely packed on the sky. Even with new techniques that 

are required to monitor flux changes of unresolved stars 
The galactic dark matter may be partly composed of com- in the face of seei ng variations JCrotte fc: Tomaney 1996| 
pact objects (e.g., faint stars, brown dwarfs, Jupiters) 'Ansari et al. 19971 lAnsari et al. 199*91 . the low signal- to- 
that reside in halos and are popularly called MACHOs noise ratio (S/N) engenders a whole range of problems. 
("MAssive Compact Halo Objects"). Microlensing sur- First, the detection efficiency is reduced. Second, there is 
veys towards M31 IjCrotts 19921 IBaillon et al. 19 931 have a degeneracy between the Einstein crossing time, the im- 
the potential to resolve the puzzling question raised by pact parameter and the source flux ( Gould 1996). Third, 
searches toward the Magellanic Clouds: the optical depth some variable stars can not be easily distinguished from 
T ~ 10~^ measured by MACHO l|Alcock et al. 2000)l is microlensing events and so will contaminate the signal, 
too large by a factor 5 to be accounted for by known pop- We elaborate on each of these points as follows: 
ulations of stars and too small by the same factor to ac- i) The loss of detection efficiency is severe because a mi- 
count for the dark matter, while the mass scale inferred crolensing event can be rejected by the selection proce- 
for the lenses M ~ 0.4 Mq is in the mid-range of normal dure if the source star or neighbouring blended stars are 
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variable. Indeed, if it is to be detected as microlensing, 
an event must rise above the photon noise due to all the 
blended neighbouring stars. For a fixed impact parameter, 
the brighter the source star, the easier it is to detect the 
event. So, bright sources are the most likely microlensing 
candidates. Unfortunately, the Hipparcos catalogue shows 
that most of the bright sources with My < are prone to 
intrinsic variability ( [Ferryman et al. 1997| ). 

ii) The degeneracy between parameters of the lightcurves 
occurs mainly around the time of maximum magnifica- 
tion and becomes more severe as the impact parameter 
becomes smaller. It can be partly broken for events with 
good S/N and good sampling on the wings - as for three 
of the four events presented later. 

iii) To distinguish between any MACHO population and 
variable stars, we intend to exploit the fact that M31 is 
highly inclined {i ^ 77°) to our line of sight. Therefore, if 
MACHOS are distributed in a roughly spherical halo, the 
density of MACHOs along the line of sight is larger on the 
far side of the M31 disk than on the near side. This im- 
plies a larger optical depth and an excess of microlensing 
events on the far side IjCrotts 19921 IKerins et al. 2001|l . 

The POINT-AGAPE collaboration is carrying out a 
pixel-lcnsing survey of M31 using the Wide Field Camera 
(WFC) on the 2.5 m Isaac Newton Telescope (INT). We 
monitor two fields, each of ^O.Sdeg^, located North 
and South of the M31 centre. After a brief description 
of the observations and data analysis in Section [3 we 
present four events with high S/N and short durations 
in Sectional for which microlensing is by far the most 
plausible interpretation. 

2. Observations and Data Analysis 

The analysed data are from 143 nights between August 
1999 and January 2001. The observations are made in 
three bands close to Sloan g',r',i'. The exposure times 
are typically between 5 and 10 minutes per night, field 
and filter. Because the total allocated time per night is 
usually less than one hour, observations are not perfomed 
in all filters each night. Moreover, the observations are 
strongly clustered in time because the WFC was not al- 
ways mounted on the telescope. 

The data reduction is described in detail by Paulin- 
Henriksson (2002) and is similar to the method given in 
previous papers IjAnsari et al. 19971 lAuriere et al. 20011 
ICalchi Novati et al. 2002|l . After bias subtraction and flat- 
fielding, each image is geometrically and photometrically 
aligned relative to reference images (one per CCD), which 
are chosen to have long exposure times, typical seeing 
between l."3 and l."6, and little contamination from the 
Moon. To remove the correlations with seeing variations, 
we first compute lightcurves on 7-pixels square "superpix- 
els". We then apply an empirical correction on the flux 
of the superpixels, called "seeing stabilisation". This is 
described briefly in Section 12.11 and will be discussed in 
more detail in a forthcoming paper. The conversion to 
Johnson/Cousins {V, R, I) is made by using the photom- 
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Fig. 1. Correlation between $f - $f and ^1°^ — ^1°^ for a 
seeing difference between image / and the reference image: 
Aseeing ^ 0.4" and (7 x 7) superpixels. 



etry standards of Haiman et al. (1994). The detection of 
events is made in the r' band, which offers the best com- 
promise between sampling and sky background. Other 
bands are then used to test the achromaticity of candi- 
dates. A bump is defined by at least three consecutive r' 
data points rising above the baseline by at least 3a. In this 
way, we detect about 80 000 variable objects. As a prelim- 
inary selection, we keep a subsample of the brightest 10%. 
More precisely, we demand R{AF) < 21, where R{AF) 
is the (Cousins) magnitude of the flux difference between 
the baseline flux and the maximally magnified flux during 
the event. Note that for small impact parameters, such as 
applies for the four candidates presented below, R(AF) is 
similar to the magnitude of the event at maximum mag- 
nification. Selection of the microlensing candidates among 
the remaining events is described in Section 12.21 

2.1. Seeing stabilisation 

For very crowded fields like ours, and in the absence of 
resolved stars, the difference between an image and its 
own median comes from star density fluctuations. This 
difference is fully correlated from image to image. The 
correlation, shown in flgure^ is 

$f-$^«(a^ + l)x ($-f-$-f) (1) 

where $f and ^l'^^ are the fluxes of superpixel i on image I 
and the reference image respectively, $i is the median flux 
computed on a 41-pixels (13."53) square centred on the 
superpixel i; and is constant over image I (zero when 
the telescope conditions are identical for image / and the 




S. Paulin-Henriksson et al.: The POINT- AG APE survey. 



3 




Aseeing (ai^csec) 

Fig. 2. The coefficient versus the seeing difference be- 
tween images and reference image for a sample of 166 im- 
ages and (7x7) superpixels. 



reference image). The slope of the correlation depends on 
the seeing difference between the two images: the larger 
the seeing on image /, the more this image is close to 
its median, and the smaller the correlation parameter . 
Figure 13 shows the correlation between and the seeing 
difference Aseeing between the image and the reference 
image. As expected, it indicates that this correlation is 
very strong. Therefore, the corrected flux 
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restores the flux of image / and superpixel i to the pho- 
tometric level of the reference image. 

Since equation |^ is a statistical correlation with an 
intrinsic width, equation ^ implies a residual gaussian 

noise which is constant over each image. The error bar on 
c I 

' is then redefined to be 



(3) 



where a^^ ^ is the photon noise in the superpixel i of im- 
age / and crfpc is the residual noise on image /. Figure 13 
shows an example of seeing stabilisation for the lightcurve 
of a typical unresolved variable star. The larger the su- 
perpixel size, the larger the noise, since a larger bin in- 
troduces more background into the data. The superpixel 
size of 7 X 7 pixels (2." 31 x 2."31) is chosen to give the 
best compromise between seeing stabilisation and S/N 
ifPaulin-Henriksson 2002,) . 

2.2. Selection of microlensing candidates 

We fit all preselected lightcurves to a standard microlens- 
ing curve fPaczyhsk i 1986) ) with seven parameters: the 
Einstein crossing time ts, the date of maximum magnifica- 
tion to, the impact parameter uq, and two flux parameters 
for each fllter, one for the source Fg and the other for the 
background Ft. To allow for non-standard microlensing 
events, we initially set a loose threshold of x^/dof < 5. 
To keep high S/N candidates, we calculate the probabil- 
ity P that a bump is due to random noise, and demand 
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Fig. 3. Lightcurve in r' filter of a 7 x 7-superpixel centred 
on an unresolved variable star before (upper panel) and 
after (lower panel) seeing stabilisation. The horizontal line 
is the baseline. 



— InP > 100 in r' and — InP > 20 in one other filter, 
with at least two points (in either band) on both the ris- 
ing and falling parts of the variation. These cuts leave 441 
candidates, one third of which show secondary peaks. Our 
purpose here is to present events for which we have very 
high confidence. To eliminate periodic variables, we reject 
by eye every lightcurve with a secondary peak comparable 
(in terms of amplitude and shape) to the microlensing can- 
didate. We then check if the remaining secondary peaks 
are due to variations in the neighbourhood with a simple 
image differencing test: for each peak, we add all images 
belonging to the peak and subtract as many images far 
from any peak and with similar seeing. On this difference 
image, variations separated by more than ^ 1" (^3 pix- 
els) are easily resolved, as shown for example in Section 
13.21 If we cannot distinguish secondary bumps from the 
microlensing candidate, the lightcurve is rejected. 

After this cull, there remain 362 candidates. To distin- 
guish microlensing events with long full width at half max- 
imum timescales (ii/2 > 25 days) from intrinsically vari- 
able stars will require additional baseline data. These will 
need to come from a third season of observations and pos- 
sibly from other telescopes, such as the MDM-McGraw- 
Hill telescope ("Calchi No vati et al. 2002'jl . Moreover, un- 
less typical MACHO masses exceed 1 Mq, we expect more 
than 80% of the microlensing events to have ii/2 < 
25 days. We therefore restrict ourselves to events shorter 
than 25 days. We expect this cut to eliminate most of the 
Mira variables. Noda et al. (2002) show that some Mira- 
like variables have large amplitudes and hence small ti/2 
and so may pass this cut. However, such variables have 
periods of the order of 100 days and so they will almost 
always have multiple bumps over our 2 years baseline. 
Then the probability for such a star to mimic a Paczyhski 
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curve is extremely small. This leaves eight candidates, four 
of which are discussed below. Of the other four events, 
one is suggestive of a binary lens and will be the sub- 
ject of further analysis. Another shows some asymmet- 
ric correlations among the residuals of the Paczyhski fit 
and is suspected to be due to a variable star rather than 
a microlensing effect. The two remaining events are not 
convincing microlensing candidates because they are too 
poorly sampled and/or have too noisy a baseline to allow 
the correlations to be studied. 



3. Four Robust Candidates 

The four candidates are PA-99-N1, PA-99-N2, PA-00-S3 
and PA-00-S4. The letter N(S) indicates whether the event 
lies in the north(south) INT WFC field, the first number 
99(00) gives the year in which the maximum occurs and 
the second number is assigned sequentially in the order 
of detection. Candidates PA-99-N1 and PA-00-S4 have al- 
ready been discussed by Auriere et al. (2001) and Paulin- 
Henriksson et al. (2002). 

The statistical relevance of the candidates is estimated 
via the total signal-to- noise ratio S /N of the bumps in the 
r' filter : 



iGbump 



(4) 



where i^bi is the baseline flux given by the Paczynski fit. 
The sum is taken over all datapoints that lie more than 
3(7 above the baseline. The four candidates have signal-to- 
noise ratios of between ~ 60 and ~ 1600, that is, high com- 
pared to typical events in the database. Tabled gives the 
main characteristics of the four candidates, while Figure 
^ shows their lightcurves. In the following, they are pre- 
sented in ascending order of distance to the centre of M31, 
as illustrated in Figure [3 



3.2. PA-99-N1 

PA-99-N1 hes 7'52" from the centre of M31 and occurred 
in August 1999, at the beginning of the first season of 
observation. The data up until November 1999 have al- 
ready been presented by Auriere et al. (2001). Since then, 
the baseline has been extended to January 2001, yielding 
better constraints on the parameters of the Paczyfiski fit 
(see Table This baseline shows two secondary bumps 
in December 1999 and November-December 2000 but, us- 
ing the simple image differencing procedure explained in 
Section O 

one can show that these secondary bumps are 
separated from the event by about 3 pixels (~ 1"), as 
shown in Figure El Therefore, the existence of these sec- 
ondary bumps is no reason to reject PA-99-N1, and data 
points belonging to secondary bumps are masked for fur- 
ther analysis. 

Our detailed study of this event (Auriere et al. 
2001) led us to conclude that the source star is almost 
certainly identified on Hubble Space Telescope (HST) 
archival images and has Johnson/Cousins magnitudes; 
/ = 22.41 ± 0.10, V = 24.51 ± 0.12. This allows one to 
break the degeneracy between the Einstein crossing time 
and the impact parameter, and so obtain direct mea- 
surements of the event duration and impact parame- 
ter: tE = 9.74 ±0.70 days, uo = 0.057 ±0.004. These val- 
ues arc slightly different from those given in Auriere et 
al. (2001), as we have subsequently extended the baseline 
for this event. If the halo fraction is above 20%, the lens 
is most probably a MACHO (with equal chance to be in 
the M31 or Milky Way halo). However, it is also plausible 
that the lens is a star, in which case the most probable 
mass is around M ~ 0.2 Mq. 



3.3. PA-99-N2 



3.1. PA-00-S3 

PA-00-S3 lies 4'00" from the centre of M31, within the 
fields observed using the 1.3 m MDM-McGraw-Hill tele- 
scope (Calchi Novati et al. 2002). The continuation of 
the lightcurve in 1998 and 1999 with data from this 
telescope shows no secondary bump, confirming the mi- 
crolensing hypothesis. At this position, bulge-bulge lens- 
ing is expected to dominate (e.g., Kerins et al. 2001). 
The high S/N of this event, together with the fact that 
the wings of the lightcurve are well sampled, has enabled 
a remarkably good direct estimate of the Einstein cross- 
ing time, ~ 13 ± 4 days. For comparison, we would ex- 
pect timescales of 20 (M/A/q)^/^ days for bulge-disk 
lensing, ~ 30 (M/Mq)^/^ days for bulge-bulge lensing or 
~ 60 {M/MqY^'^ days for disk-disk lensing. The proximity 
of PA-00-S3 to the bulge, together with its Einstein cross- 
ing time, strongly suggests that it is a low mass stellar 
lens event. 



At 22'03" from the centre of M31, PA-99-N2 hes weU 
outside the projected area of the M31 bulge and there- 
fore appears to be an excellent candidate for a MACHO 
lens. Precise determination of all five parameters of the 
Paczynski fit is possible because of very good sampling 
and small error bars. We find a duration for this events of 
tE 92 ± 4 days and maximum magnification ^max = 
13.3 ± 0.7. MACHOs are expected to have a duration 
centred on ~ 50 (M/M©)^/^ days. The alternative to 
MACHOs is self-lensing by an M31 disk star l|Gould 1994|l . 
with a typical duration Ie 60 (M/M©)^/^ days. 

The very high S/N allows achromatic deviations to 
the standard Paczyhski curve to be revealed, implying a 

per degree of freedom of 3.1. Some evidence of this is 
visible in Figure 0] If these deviations are not simply due 
to systematic photometric errors, they are suggestive of a 
parallax effect (Alcock ct a l. 1997|l and/or a close caustic 
approach ((Albrow et al. 2002)1 . 
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Fig. 4. Lightcurves of the four microlensing candidates (presented in Sections 1^.11 to rO|l . For each event, the top panel 
shows two seasons of analysed data in the r' filter. Lower panels are zooms that focus on candidates in all bands for 
which data have been taken. Solid lines are best-fit Paczyhski (1986) curves. For PA-99-N1, secondary bumps due to a 
neighbouring variable star are masked (see Section IX^ . Note that deviations from the Paczyhski curve for PA-99-N2 
are achromatic (see Section rO)l . 



3.4. PA-00-S4 

We have discussed previously PA-00-S4 (Paulin- 
Henriksson et al. 2002). It lies 22'31" from the centre of 
M31, but unlike PA-99-N2, there is a concentration of 
stellar lenses along the line of sight. Indeed, it is only 2'54" 
from the centre of M32. The blue color {R-I ^ 0.0±0.14) 
of this event argues strongly for a source that lies in 
the disk of M31. Because M32 is believed to be in the 



foreground ( |Byrd 1976| |Ford, Jacoby and Jenner 1978| ), 
the proximity of the lens to the line of sight of M32 
suggests that it is most naturally interpreted as a star 
in M32. However, this argument is not conclusive and a 
MACHO interpretation is still possible. 
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PA-99-N1 


PA-99-N2 


PA-00-S3 


PA-00-S4 


Q (J2000) 


00h42m5 1.42s 


00M4m20.81s 


00h42m30.51s 


00h42m29.97s 


S (J2000) 


41°23'53.7" 


41°28'45.2" 


41° 13' 04.9" 


40°53'47.1" 


d 


7' 52" 


22'03" 


4'00" 


22'31" 


ii/2 (days) 


1.80 ±0.22 


21.75 ± 0.20 


2.18 ±0.14 


2.09 ±0.11 


to (JD-2451392.5) 


13.87 ±0.06 


71.56 ±0.10 


458.35 ± 0.02 


488.90 ± 0.07 




0.057 ± 0.004 


0.075 ± 0.004 


0.053lH?6 


0.00472+»og«J« 




17.54ti;?^ 


13.33+!]-^? 


18.88l|;8| 


211+i^r 


tE (days) 


9.74 ± 0.70 


91.911^83 


1 2 . 56 2 23 


128.581^^2^?! 


AR 


20.80 ±0.13 


19.0 ± 0.2 


18.8 ± 0.2 


20.7 ±0.2 


V ~R 


1.2 ±0.2 


1.0 ±0.1 






R~I 






0.6 ±0.1 


0.0 ±0.1 


S/N 
xVd.f. 


63.6 (7 pts) 


1603.0 (54 pts) 


115.7 (5 pts) 


116.2 (10 pts) 


0.9 (65 d.f.) 


3.1 (233 d.f.) 


1.1 (172 d.f.) 


0.7 (146 d.f.) 



Table 1. Main characteristics of the four microlcnsing candidates, d: distance from the centre of M31. ti/2'- full width 
at half maximum timescale. date of maximum magnification, uq^ impact parameter in units of the Einstein radius, 
^max^ maximum magnification I/mq for such large magnifications), Einstein crossing time. For PA-99-N2, PA- 
00-S3 and PA-00-s4, mq, A„iax and t-E, are computed without any prior on the source fluxes. The source star of PA-99-N1 
has been identified on HST archival images (Auriere et al. 2001). Ai?: Cousins magnitude of the maximum source 
flux variation, which is also the magnitude at tg for such large magnifications. V", R, I: Johnson/Cousins magnitudes. 
S/N: total signal-to-noise ratio, as defined in equation^ the number quoted in parenthesis indicates the number of 
points (in r' filter) contributing to S/N. x^/d.f.: per degree of freedom of the Paczyhski fit; the number quoted in 
parenthesis indicates the number of degree of freedom of the fit. 



4. Conclusion 

We have reported first results extracted from two years of 
data in our pixel lensing survey of the Andromeda galaxy. 
By imposing the stringent requirements i?(Ai^) < 21 and 
ti/2 < 25 days, we have selected four very convincing mi- 
crolensing candidates with high signal-to-noise ratio and 
short timescales. For three of our four events we have been 
able to make reliable determinations of the Einstein cross- 
ing time, which provides additional clues as to the prob- 
able origin of these events. In the case of PA-00-S3, the 
event is most likely caused by a stellar lens in the M31 
bulge. In the three other cases, MACHOs and stellar lens- 
ing are plausible. 
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